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We studied the response of the Josephson vortex lattice in layered superconductors to the high- 
frequency c-axis electric field. Wo found a simple relation connecting the dynamic dielectric constant 
with the perturbation of the superconducting phase, induced by oscillating electric field. Numerically 
solving equations for the oscillating phases, we computed the frequency dependences of the loss 
function at different magnetic fields, including regions of both dilute and dense Josephson vortex 
lattices. The overall behavior is mainly determined by the c-axis and in-plane dissipation parameters, 
which are inversely proportional to the anisotropy. The cases of weak and strong dissipations are 
realized in Bi2Sr2CaCu2 0x and underdoped YBa2Cu30x, respectively. The main feature of the 
response is the Josephson-plasma-rcsonancc peak. In the weak-dissipation case additional satellites 
appear in the dilute regime in the higher-frequency region due to the excitation of the plasma 
modes with the wave vectors set by the lattice structure. In the dense-lattice limit the plasma 
peak moves to a higher frequency and its intensity rapidly decreases, in agreement with experiment 
and analytical theory. The behavior of the loss function at low frequencies is well described by 
the phcnomcnological theory of vortex oscillations. In the case of very strong in-plane dissipation 
an additional peak in the loss function appears below the plasma frequency. Such peak has been 
observed experimentally in underdoped YBa2Cu3 0x. It is caused by the frequency dependence of 
the in-plane contribution to losses rather than a definite mode of phase oscillations. 



I. INTRODUCTION 

The Josephson plasma resonance ( JPR) is one of the 
most prominent manifestations of the intrinsic Josephson 
effect in layered superconductors.^'^ It has been estab- 
lished valuable tool to study the intrinsic proper- 
tics of superconductors®'^ and it was extensively used to 
probe different states of vortex matter. The value of 
the JPR frequency depends on the anisotropy factor and 
it ranges widely for different compounds, from several 
hundred gigahertz in Bi2Sr2CaCu20x (BSCCO) to sev- 
eral terahertz in underdoped YBa2Cu30x (YBCO). 

An interesting issue is the influence of the magnetic 
field applied along the layer direction on the high- 
frequency response and, in particular, on the JPR. Such 
field forms the lattice of Josephson vortices (JVL). The 
anisotropy factor, 7, and the interlayer periodicity, s, set 
the important field scale, Bcr = $0/(27775^). At B < Ba 
the Josephson vortices are well separated and form a di- 
lute lattice. When the magnetic field exceeds B^i the 
Josephson vortices homogeneously fill all layers (dense- 
lattice regime). The crossover field ranges from ~0.5 
tesla for BSCCO to ~10 tesla for underdoped YBCO. 
The JVL state is characterized by a very rich spectrum of 
dynamic properties. In particular, the transport proper- 
ties of the JVL in BSCCO have been extensively studied 
by several experimental groups. ^^"^"^ 

The high-frequency response in the magnetic fields 
along the layers has been studied experimentally in 
BSCCO using the microwave absorption in cavity 
resonators^^, in La2_xSrxCu04^, and in underdoped 
YBCO^^'^'^ by the infrared reflection spectroscopy. A 
detailed comparison between the behaviors of the high- 
frequency response in the in-plane magnetic field for 
these compounds has been made recently by LaForge 



et al?"^ In BSCCO, two resonance absorption peaks 
have been observed. ^-'^ The upper-resonance frequency in- 
creases with the field and approaches the JPR frequency 
at small fields, while the lower-resonance frequency de- 
creases with the field and approaches approximately half 
of the JPR frequency at small fields. In the underdoped 
YBCO^^'^^, the JPR peak in the loss function rapidly 
loses its intensity with increasing the field, while the res- 
onance frequency either does not move or slowly increases 
with the field. In addition, another wide peak emerges 
at a smaller frequency and its intensity increases with 
increasing field. 

Several theoretical approaches have been used to de- 
scribe the response of the JVL state to the oscillating 
electric field. A phenomenological vortex-oscillation the- 
ory has been proposed by Tachiki et alp- . This theory de- 
scribes the response of the JVL at small frequencies and 
fields, in terms of phenomenological vortex parameters, 
mass, viscosity coefficient, and pinning spring constant. 
This approach is expected to work at frequencies much 
smaller then the JPR frequency, i.e., it can not be used 
to describe the JPR peak itself. 

On the other hand, the plasmon spectrum at high mag- 
netic fields, in the dense-lattice limit, has been calculated 

by Bulaevskii et al?'^ It was found that the plasma mode 
at zero wave vector increases proportionally to the mag- 
netic field, ujp{B) = ujp{0)B/{2Bci)- This prediction de- 
scribes very well the behavior of the high-frequency mode 
in BSCCO. A more quantitative numerical study of the 
JVL plasma modes in the dense-lattice regime has been 
done by Koyama.^^ He reproduced the mode, linearly 
growing with field, and also found additional modes at 
smaller frequencies. 

In this paper we develop a quantitative description of 
the high-frequency response for a homogenous layered 
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superconductor valid for whole range of frequencies and 
fields. We relate the dynamic dielectric constant via sim- 
ple averages of the oscillating phases. The high-frequency 
response is mainly determined by the c-axis and in-plane 
dissipation parameters, which are inversely proportional 
to the anisotropy. Analytical results for the dynamic 
dielectric constant and loss function can be derived in 
limiting cases, at small fields and frequencies (vortex- 
oscillation regime) and at high fields, in the dense- lattice 
regime. Numerically solving dynamic equations for the 
oscillating phases, we studied applicability limits of the 
approximate analytical results and investigated the in- 
fluence of the dissipation parameters on the shape of the 
loss function. Computing the oscillating phases for dif- 
ferent vortex lattices, we traced the field evolution of the 
loss function with increasing magnetic field for the cases 
of weak and strong dissipations, typical for BSCCO and 
underdoped YBCO, respectively. 



II. DYNAMIC PHASE EQUATIONS, 
DIELECTRIC CONSTANT, AND LOSS 
FUNCTION 

The equations describing phase dynamics in layered su- 
perconductors can be derived from Maxwell's equations 
expressing fields and currents in terms of the gauge in- 
variant phase difference between the layers On = 4>n+i — 
(jin — (27rs/$o) ^z- These equations have been presented 
in several different forms. We will use them in the form 
of coupled equations for the phase differences and mag- 
netic field when charging effects can be neglected, see, 
e.g., Refs. 27,28, 
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where the magnetic field is along the y axis, aab and 
CTc are the components of the quasiparticle conductivity, 
Aab and Ac are the components of the London penetra- 
tion depth, jj = c'l'o/(87r^,sA^) is the Josephson current 
density, Up = cj y^Ac is the plasma frequency, is the 
external electric field, and V^Bn = -Bn-i-i + -Bn-i — 2B„. 
Neglecting charging effects^'', the local electric field is 
connected with the phase difference by the Josephson re- 
lation, 
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The average magnetic induction inside the superconduc- 
tor, By^ fixes the average phase gradient, (d9n ^/dx^ = 
2-KsBy/^o. 

To facilitate analysis, we use a standard transformation 

to the reduced variables, 



with Aj = 7S, and introduce the dimensionless parame- 
ters 
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I ■ is important to note that both damping parameters, 
Vc and Vab, roughly scale inversely proportional to the 
anisotropy factor 7, meaning that the effective damping 
is stronger in less anisotropic materials. In particular, 
the cases of weak and strong dissipation are realized in 
BSCCO and underdoped YBCO, respectively. Due to the 
d-wave pairing in the high-temperature superconductors, 
both dissipation parameters do not vanish at T ^ 0. 
Another important feature of the high-temperature su- 
perconductors is that the in-plane dissipation is typically 
much stronger than the c-axis dissipation, Vab ^ ^c-^^ 
This is a consequence of a rapid decrease of the in-plane 
scattering rate with decreasing temperature, which mani- 
fests itself as a large peak in the temperature dependence 
of the in-plane quasiparticle conductivity. ^^'^^ 

Assuming an oscillating external field and using a com- 
plex presentation, 'Dz{t) = exp(— iwt), we obtain for 
small oscillations 
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where & ~ t^/wp, C„(x) = cos[6n\x)], and the static 

phases, 9n \x), are determined by the following reduced 
equations 
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with the addition condition {^dOn^ /dx^ — h = 
27r7S^i?y/$o- We introduce the reduced oscillating phase 
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for which we can derive from Eqs. (3a) and (3b) the fol- 
lowing reduced equation 
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Prom the Josephson relation (2) we find E^, = 
{—iui/Ec) '&T^z, meaning that the dynamic dielectric con- 
stant is connected with the average reduced oscillating 
phase, 1?, by a simple relation 
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Consider the case of a zero magnetic field first. In this 
case On^ {x) = and the oscillating phase is given by 



—uv 



(7) 



In this case Eq. (6) gives well-known results for the dy- 
namic dielectric constant, eco(w) = Vz/Ez, and the loss 
function, Lq{u>) = Im[— l/£co(w)] at the zero magnetic 
field, which we present in real units. 
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The zero-field loss function has a peak at the Josephson 
plasma frequency with width determined only by the c- 
axis quasiparticle conductivity. 

Consider now the case of a finite magnetic field applied 
in the layer direction. Such magnetic field generates the 
lattice of the Josephson vortices, which is described by 
the static phase distribution, On\x), obeying Eq. (4). 
The oscillating phase, in turn, is fully determined by this 
phase via the spatial distribution of the average cosine, 
Cn{x). Averaging Eq. (4), we obtain an obvious identity 

(sin^i^^) = 0, indicating that there is no average current 
in the ground state. In the field range By S> $o/27rAa(,Ac 
the term 1//^ can be neglected. In this limit it is more 
convenient to operate with the in-plane phases (f>n\x)^ 



defined by the relation 6n^ - 
obey the following equation 
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Due to the translational invariance, every solution (j)n^ (x) 

generates a continuous family of solutions 4>li\x — u) 
corresponding to the lattice shifts u. In particular, the 
phase change for a small displacement is given by = 
/dx. Taking the a;-derivative of Eq. (10), we find 
that d(pn^ /dx obeys the following equation 
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with Cn = cos (4>n+l ~ 4>n^ + ^"^^ V„C„ = C„ - 

Cn-1- The condition {^d sin 9n^ /dx^ — gives another 
useful identity for d4in^ /dx 
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We now proceed with the analysis of the dynamic 
phase equation (5). Averaging this equation, we obtain 
the following relation 



[iUcU + W^) 1? — Cn{x)'dn = il^- 

Splitting ■!?„(x) into the average and oscillating-in-space 
parts, 'i?„(a;) ='&{! + w„(a;)), with Wn{x) = 0, we obtain 
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Using relation (6), we obtain the following result for the 
dielectric constant 



a;2 



ec('^) Cb'^ - C - Cn{x)wn + iyci^ 



(13) 



Therefore, the dynamic dielectric constant is fully deter- 
mined by the two simple averaged quantities, the static 
average cosine, C, and the average i ncluding the spatial 
distribution of the oscillating phase, C„(a:)w„. 

Introducing again the in-plane oscillating phases, u>„ = 
4>n+i — 4>ni and neglecting terms of the order of l/P-, we 
derive the following equation 
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The oscillating phase has the same symmetry properties 
as the vortex lattice. Comparing this equation with Eq. 
(11), we see that in the limit w ^ the solution is given 

by 



id_^ 

h dx 



(15) 



This solution corresponds to the homogeneous lattice 
shift. As follows from Eq. (12), the combination C + 

Cn{x)Wn = C + Cn{x) {4>n+i " 4'n)i which determines 
the dielectric constant (13), vanishes in this limit. This 

property is a consequence of the translational invariance 
and it is only true in the absence of pinning of the vortex 
lattice. 

In summary, to find the dynamic dielectric constant at 
a given field, one has to find first the static phase from Eq. 
(10) assuming a definite vortex-lattice structure and com- 
pute the average cosine, C = (^cos (^4>n+i ~ ^rf"* + hx^ ^. 
After that, one has to solve the dynamic equation (14) 
and compute the average C„(a:)«;„ = C„(x) (ci„+i— 0„). 
These two averages completely determine £, (0;) via 
Eq. (13). In the following sections, we will consider 
regimes in which analytical solutions are possible, the 
high-field limit and the vortex-oscillation regime at small 
frequencies. Then, we will present the results of the nu- 
merical analysis in the full field and frequency range in 
the cases of weak dissipation (BSCCO) and large dissi- 
pation (underdoped YBCO). 
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III. HIGH-FIELD REGIME 

In this section we consider the high-field regime, ft, ^ 1, 
corresponding to the dense-lattice limit in which all in- 
terlayer junctions are homogeneously filled with vortices. 
This regime allows for the full analytical description using 
an expansion with respect to the Josephson currents. In 
particular, the spectrum of the plasma modes and their 
damping in this limit have been found by Bulaevskii et 
al.'^^ We use the same approach to derive the dynamic 
dielectric constant. 

The static phase solution at high fields, describing the 
triangular lattice, is given by 



7rn(n — 1) 



■ sin {hu + TTn) . 
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Using this distribution, we compute the average cosine 
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At high fields we can neglect rapidly oscillating C„(a;) in 
the right hand side of Eq. (14). This allows us to obtain 
the solution 
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compute the average 

Cn{x) (0n+l - 0,1) 
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and, finally, obtain the high-field limit of E(\lo) from Eq. 
(13), 
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In the low-damping regime v^^ v^h ^ 1 the loss function 
has a peak at cj = /i/2 (in real units uj — ujpiT^s^By/^t)) 
corresponding to the homogeneous plasma mode.^"* Such 
linear growth of the plasma frequency with field has been 
indeed observed experimentally in underdoped BSCCO 
by Kakeya et al.^^ 

To verify the accuracy of the high-field approximation 
(19), we compare in Fig. 1 the loss function obtained 
from this formula with an accurate numerical solution 
for two field values, ft, = 2 and 4 in the case of weak 
dissipation, Vc — 0.01 and i^ab = 0.1. One can see that the 
high-field formula accurately describes the low-frequency 
region already at ft = 2, but it overestimates the peak 
frequency. At ft = 4 the high-frequency approximation is 
already undistinguishable from the exact result. 



IV. VORTEX-OSCILLATIONS REGIME AT 
SMALL FIELDS AND FREQUENCIES 

In this section we consider the phenomenological the- 
ory of vortex oscillations which describes the response 
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FIG. 1: (Color online) Comparison between the approximate 
high-field limit of the loss function based on Eq. (19) and 
exact numerical solution for two field values, ft, = 2 and 4 in 
the case of weak dissipation, — 0.01 and i/at ~ 0.1. For 
clarity, the vertical scale for the ft = 4 plots is magnified 5 
times. 



of the vortex lattice at small frequencies, lo cOp. For 
the Abrikosov vortex lattice such a theory was developed 
by Coffey and Clem.'^° The dynamic dielectric constant 
for the Josephson vortex lattice has been derived using a 
similar approach by Tachiki, Koyama, and Takahashi.^ 

Consider a superconductor in the vortex state carrying 
ac supercurrent js oc exp(— iwi) along the c axis. The ac 
electric field consists of the London term and the contri- 
bution from the vortex oscillations. 
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The vortex oscillation amplitude, 
the equation 



can be found from 



-pjuj - irijuj - 



K) u 
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where pj is the linear mass of the Josephson vortex^**, 
rjj is its viscosity coefficient'^^'^^, and K is the spring 
constant due to pinning (it is neglected in the rest of the 
paper). In Appendix A we present formulas for pj and 
rjj in terms of superconductor parameters. Finding the 
oscillating amplitude, we obtain 



A-kK 



—pjuj^ — irjjuj + 



K 
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A finite electric field also generates the quasiparticle cur- 
rent jn — (Tc.7i{'^)Ez- Therefore, the total conductivity 
ac{By,uj) = {is + in)/E:, is given by 



ac{By,uj) =crc.n(w) 
,2 



p 



1- 



1 



By<^0 



—pjuj^ — irjjuj+ 



K AttXI 



(23) 



5 



The interesting feature of the real part of conductivity, 
ac,i{By,u}) = Re[(Tc,„(w)] 

rijEcUjlBy^ol (47rAc)^ 



Role of in-plane dissipation 



Role of c-axis dissipation 



(24) 

is a resonance peak at the frequency ujr — 
y^[K + By^Q/ (AttXI)] Ipj. As one can see from Eq. (22), 
this resonance is a result of compensation of the Lon- 
don and vortex contributions to the oscillating electric 
field at Using the result for the vortex linear 

mass pj from Appendix A and neglecting pinning, the 
formula for lo^ can be rewritten more transparently as 
ujr ~ 0.84cjp^27r7s2i3y/<i)o. At w 0, Eq. (24) gives a 
known result for the dc flux-flow conductivity."^^ 

Using the well-known relation between the dynamic 
dielectric constant and conductivity, ec(w) = £c ~ 
47rcrc(w)/iu;, we obtain^ 



EciBy^Lo) =ec 



47r(Tc,n(w) 



By<i>0 
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Using again formulas for pj and 77,7 from Appendix A 
and neglecting pinning, we rewrite this result in reduced 
variables in the form, convenient for comparison with nu- 
merical calculations, 



£c{h,Lo) 



1^ 



1/w 



~,2 



UJ l — 2'Kh/[Cc{oj^+il^c>Jj)+CabiVabOj] 

(26) 

with Cc ~ 9.0 and Cab ~ 2.4. We remind that this for- 
mula is valid only at small fields, /i <C 1, and frequencies, 
oX 1. 



V. FIELD EVOLUTION OF THE LOSS 
FUNCTION 

In the full range of frequencies and fields, the dynamic 
dielectric function can only be computed numerically. At 
the first step, one has to find the static phase distri- 
bution from Eq. (10) assuming a definite vortex lattice 
structure. To probe general trends, we limit ourselves 
here only by simple aligned lattices. At a fixed mag- 
netic field, such a lattice is fully defined by the num- 
ber of layers separating the layers filled with vortices, 
Nz, see sketch of the aligned lattice with A^^ = 2 in 
the inset of Fig. 3. At small fields, such lattices are 
realized in ground states in the vicinity of two sets of 
commensurate fields, Bi{Nz) = V^^o/i'^N^ls'^) and 
B2{N,) = <^>o/{2V^N^js^) {hi{N,) = ttV3/N^ and 
h2{Nz) — 7r/(\/3-/V^) in reduced units). At the interme- 
diate field values the ground state is given by misaligned 
lattices, i^'i"^ 




FIG. 2: (Color online) Influence of the dissipation parame- 
ters on the loss-function shape in the dilute-lattice regime at 
h = 0.5 and = 2. In-plane dissipation (left plot) does 
not influence much the JPR peak, but strongly influences the 
shape below the peak. 



At the first stage, we solved the static phase equa- 
tions (10) for fixed h and N^. This solution has been 
used as an input for the dynamic phase equations (14). 
Finally, the oscillating phase determines the dynamic 
dielectric constant (13) and the reduced loss function 
L{u}) = -Im[ec/£c('^)]- 

We start with the discussion of several general prop- 
erties of the high-frequency response. To illustrate the 
roles of two different dissipation channels, we show in 
Fig. 2 series of the frequency dependencies of the loss 
function when only one damping constant is changed for 
fixed field h — 0.5 and ~ 2 corresponding to a dilute 
lattice. We can see that the width of the JPR peak is 
mainly determined by the c-axis dissipation, while the 
in-plane dissipation has no visible infiuence on the loss 
function near the JPR peak. On the other hand, the in- 
plane dissipation strongly infiuences the shape of the loss 
function below the peak. In particular, at a very strong 
in-plane dissipation a peaklike feature appears below the 
JPR frequency which looks like an oscillation mode. 

In Fig. 3 we compare the behavior of the numerically 
computed loss function at low frequencies with predic- 
tions of the vortex-oscillation model described in section 
IV. We can see that this model typically accurately de- 
scribes the behavior at frequencies smaller then half the 
plasma frequency. At low dissipation, the loss function 
has an additional peak at small frequencies. This peak 
disappears with increasing dissipation. We have to em- 
phasize that these plots are made for an ideal homoge- 
neous system. In a real layered superconductor, pinning 
of the Josephson vortices by inhomogeneities will lead to 
the resonance peak at the pinning frequency. 
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FIG. 3: (Color online) Comparison between the vortex- 
oscillation model and exact numerical calculation for the loss 
function at several representative values of the dissipation pa- 
rameters, field, and N^, (shown in the plots). One can see that 
this model typically describes the high-frequency response 
roughly up to half of the plasma frequency. The inset in 
the upper plot illustrates the aligned lattice with A^^ = 2. 



A. Small dissipation (BSCCO) 

To illustrate a theoretically expected behavior of the 
high-frequency response in BSCCO, we made calcula- 
tions with small dissipation parameters, typical for this 
compound, Vc — 0.01 and Vab — 0.1. Figure 4 shows the 
evolution of the loss function with increasing magnetic 
field in the frequency range near the plasma resonance. 
For such small dissipation, the loss function at zero field 
has a very sharp peak at the plasma frequency. In the 
dilute-lattice regime at small fields, this peak decreases in 
amplitude and is displaced to slightly lower frequencies. 

The most interesting feature of the dilute regime is the 
appearance of the satellite peaks in the high-frequency 
part. The strongest satellite is observed at A'^ = 2 near 
h ~ 0.45. A physical origin of these satellites is clear. In 
the vortex lattice state, the homogeneous oscillations are 
coupled to the plasma modes with the wave vectors set 
by the lattice. Therefore, the location of the satellites is 
determined by the plasma spectrum. For reference, we 
present the plasma frequencies at the reciprocal-lattice 
vectors in Appendix B. At small dissipations the shape 



of the loss function is not rigidly fixed by the value of 
the magnetic field, it is also sensitive to the lattice struc- 
ture. This is illustrated in Fig. 5, where the loss function 
is plotted at fixed h for two values of N^, for which the 
lattice energies are close. One can see that the loss func- 
tion is almost Nz independent at frequencies below the 
JPR peak indicating that in this frequency range the vor- 
tices contribute independently to the response. However, 
above the peak the shape of the loss function is very sensi- 
tive to Nz ■ This means that the high-frequency response 
can potentially be used to probe the lattice structure. 

In the dense-lattice limit, the plasma peak moves to a 
higher frequency and its intensity rapidly decreases, as 
predicted by the analytical theory in Section III. This 
behavior is also consistent with experiments^. The tran- 
sition to the dense lattice is very distinct. While for 
Nz = 2 and h — 1 the peak is still located very close to 
the zero-field JPR frequency and its amplitude is compa- 
rable with the zero-field peak, at very close field h = 1.4 
in the dense-lattice regime, A^^ = 1, the peak is notice- 
ably shifted to a higher frequency, its amplitude consid- 
erably dropped, and the width increased. 

An additional broad peak exists at low frequencies. Its 
intensity is much smaller than the JPR peak. The evolu- 
tion of this peak with increasing magnetic field is shown 
in the inset of Fig. 4. One can see that the intensity 
of this peak monotonically increases with the magnetic 
field. This peak is well described by the vortex-oscillation 
model. 

We did not find any intrinsic resonances in the loss 
functions at frequencies smaller than the JPR frequency, 
meaning that there are no modes coupled to homoge- 
neous oscillations in this frequency range for a homoge- 
neous superconductor."^^ This suggests that the resonance 
feature observed in underdoped BSCCO at ~ 0.5wp in 
Ref. 21 probably has an extrinsic origin, e.g., it may be 
caused by the pinning of the Josephson vortices. 



B. Large dissipation (underdoped YBCO) 

In this subsection we discuss the high-frequency re- 
sponse of the Josephson vortex lattice in the case of large 
effective dissipation, which is realized, e.g., for the under- 
doped YBCO.SS'S"^ Figure 6 shows the frequency depen- 
dences of the loss functions for the dissipation parameters 
Vc = 0.32 and Vab — 6.0. For comparison, we also present 
the loss functions at different in-plane fields for an un- 
derdoped YBCO sample from Ref. 22. One can see that 
two experimental features are well reproduced by the the- 
ory: decreasing intensity of the main peak with increas- 
ing field and the appearance of the satellite peak in the 
low-frequency part of the line. The origin of this satel- 
lite peak is distinctly different from the high-frequency 
satellites found for small dissipation. Figure 2 demon- 
strates that this peak only appears for sufficiently strong 
in-plane dissipations. As this peak is absent in the low- 
dissipation limit, it does not correspond to any specific 
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FIG. 4: (Color online) Series of the frequency-dependent loss functions near JPR frequency at different fields for — 0.01 and 
i^ab = 0.1. The inset shows the low-frequency region for the same set of parameters. 




FIG. 5: (Color online) Plots of the loss function at the same 
h but with different for i/c = 0.01 and Vab =0.1. This 
figure illustrates the sensitivity of the loss-function shape to 
the vortex lattice structure (value of A''^) above the JPR peak. 



mode of phase oscillations. It is caused by the decrease 
of the relative contribution of the in-plane dissipation 
channel with increasing frequency. To verify this inter- 
pretation, we plot in Fig. 7 the loss function for ~ 2 
and h — 1 together with the relative contribution of in- 
plane dissipation to losses, Rab — o'abE^/iccE^ + crabE^), 
which can be rewritten in reduced coordinates. 



Rab — 



yab\d(t)n/dx\' 



1 + Wn+1 



Vab\dcj)n/dx\ 



(27) 



One can see that Rab decreases with increasing frequency 
and almost vanishes above the plasma resonance. The 
small-frequency peak roughly corresponds to the fre- 
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FIG. 6: (Color online) (a)Series of the frequency-dependent 
loss functions at different fields for = 0.32 and Vab = 6.0. 
For comparison plot (b) shows the evolution of the loss func- 
tion with increasing in-plane field for underdoped YBCO from 
Ref. 22. 



quency where Rab drops to one half. The decrease of 
the in-plane dissipation also leads to the peak in the real 
part of the optical conductivity cri, also shown in Fig. 7. 
Note that the peak in cti corresponds to the dip in the 
loss function. 

At high fields, in the dense-lattice limit, only a sin- 
gle broad dissipation peak remains. In contrast to the 
low-dissipation case, this peak is displaced to lower fre- 
quencies with increasing magnetic field. This behavior is 
not verified experimentally yet, data are only available 
for field range h < 1.22,23 
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v,=0.32, V3b=6.0, N,=2, h=1 




dependence of the in-plane contribution to losses. 
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APPENDIX A: JOSEPHSON- VORTEX MASS 
AND VISCOSITY 

The viscosity of an isolated Josephson vortex has been 
considered in Refs. 35 and 36. In the case when the dis- 
sipation is caused both c-axis and in-plane quasiparticle 
transport, the viscosity coefficient is given by 



Vj 



7S^ 



^0 

2nc 



C„ 



' 1 
1 



2 \CcVc 

TT (47rcs) 7 



CahVah\ 



(Al) 



FIG. 7: (Color online) The frequency dependences of the loss 
function, real part of conductivity, and relative contribution 
of in-plane dissipation (27) for the same damping parameters 
as in Fig. 6, N^, = 2, and h — \. Rapid decreasing of in-plane 
inhomogeneity of the oscillations with increasing frequency 
can also be seen from the visualization of oscillating patterns 
of the local electric fields at different frequencies.^* 



VI. SUMMARY 

In summary, we developed a comprehensive description 
of the high-frequency response of the Josephson vortex 
lattice in layered superconductors. We found the gen- 
eral relation (13) connecting the dynamic dielectric con- 
stant with the averages containing the static and oscil- 
lating phases. Analytical formulas for the dynamic di- 
electric constant and loss function were derived for the 
high-field regime and the vortex-oscillation regime at low 
frequencies. Numerically solving equations for the os- 
cillating phases, we explored the evolution of the loss 
function with increasing magnetic field for the cases of 
weak dissipation describing BSCCO and strong dissipa- 
tion describing underdoped YBCO. Several new features 
were found. The most interesting feature in the weak- 
dissipation case is the high-frequency satellites in the 
dilute-lattice regime, which appear due to the excitation 
of plasma modes at the wave vectors of the reciprocal lat- 
tice. In the strong-dissipation limit, we reproduced the 
additional peak in the loss function below the JPR peak 
experimentally observed in underdoped YBCO.^^'^'^ We 
established that this peak appears due to the frequency 



where the numerical constants Cc and Cah are determined 
by the phase distribution of an isolated Josephson vortex. 




The mass of the Josephson vortex has been considered 
in Ref. 34. This mass is determined by the kinetic energy, 
which for a moving Josephson vortex can be presented 

as 




For a slowly moving vortex with velocity u, the phase dif- 
ference is determined by its static distribution, &„(r, t) = 

d'^rt' (x — vt) , and we obtain 

As, by definition, £k = Lypjv'^ /2, we obtain the following 
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result for the linear vortex mass 



Pj 



7r7 [AttcsY 



oo 

E 

n=— oo 



dx 



dx 



(A2) 



The reduced combination in the dynamic dielectric con- 
stant (25) can be represented as 



2nh 



APPENDIX B: PLASMA FREQUENCIES AT 
RECIPROCAL LATTICE VECTORS: EXPECTED 
LOCATION OF SATELLITES 

The satellite peaks are expected approximately at the 
plasma frequencies for the wave vectors of the reciprocal 

lattice. In the reduced coordinates, the plasma spectrum 
at the zero magnetic field and without the charging ef- 
fects is given by 

k 



where k and q arc the wave- vector components along and 
perpendicular to the layers respectively. The reciprocal 
vectors of the vortex lattice are given by 



{k{n,m),qin,m)) = (27rn/a, 27r(m - n/2)/N^), 



where n and m arc the integer indices and a is the vortex 
separation, a = 2Tr /hN^- Neglecting the small factor 
l/l^, we obtain for the satellite frequencies 



' 1 - cos 27r (m - n 2) ^ ' 



The most intense satellite is expected for the basic recip- 
rocal vector with indices (n, m) = (1,0) 



'^'(1,0) = 1 



{hN,f 12 
1 — cos [k/Nz 



w^(k, q) = l+ iip^2{l-cosq) ^^^^ particular, for N, = 2, this gives a;2(l, 0) = 1 + 2h? 
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